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Sumraary/Overview 

Experimental,  theoretical,  and  computational  studies  of  various  combustion  phenomena  were 
carried  out  for  a  wide  range  of  conditions.  The  parameters  that  were  considered  include  the  fuel 
type,  reactant  composition,  flame  temperature,  pressure,  and  combustion  mode.  Flames  of 
hydrogen,  carbon  monoxide,  as  well  as  neat  and  mixtures  of  gaseous  and  liquid  hydrocarbons 
were  studied  experimentally  in  the  counterflow  configuration.  The  experimental  data  were  used 
as  the  basis  for  developing  and  validating  reaction  and  molecular  transport  models. 

During  the  reporting  period,  the  following  tasks  were  completed  (numbers  correspond  to  the 
archival  publications  that  are  shown  in  Section  7.0  of  Technical  Discussion): 

1.  Development  of  a  new  non-linear  extrapolation  technique  for  the  determination  of 
laminar  flame  speeds  [I]; 

2.  Molecular  transport  and  kinetics  effects  on  flame  propagation  and  extinction  in  liquid 
hydrocarbon  flames  [2]; 

3.  Quantification  of  the  effect  of  fuel  blending  on  the  flame  response  [3,4]; 

4.  Derivation  of  an  extensive  experimental  flame  database  for  Cs-Cn  n-alkane, 
cyclohexane,  and  mono-alkylated  cyclohexane  flames  [1,5]; 

5.  Fundamental  kinetic  studies  of  elementary  reactions  critical  to  a  predictive  combustion 
reaction  model  [6,7]; 
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6.  Refinement  of  USC  Mech  II  to  predict  H2/CO  combustion  at  pressures  up  to  500  atm  [8]; 

7.  Development  of  detailed  and  simplified  models  for  the  oxidation  of  Cs-Cu  «-alkanes 
[1,9,10]; 

8.  Development  of  the  method  of  uncertainty  minimization  by  polynomial  chaos  expansion 
(MUM-PCE)[lll; 

9.  Gas-kinetic  theory  analysis  of  scattering  between  gas  and  molecular  clusters  with 
application  in  the  transport  properties  of  high  molecular  weight  hydrocarbons  in  dilute 
gases  [12]. 


Technical  Discussion 
1.0  Background 

The  three-year  research  program  that  was  completed,  aimed  to  provide  insight  into  the 
physicochemical  processes  that  control  the  burning  behavior  of  fuels  that  are  of  relevance  to 
high-speed  air-breathing  propulsion.  This  was  achieved  through  combined  experimental, 
theoretical,  and  detailed  numerical  studies  of  a  number  of  fundamental  combustion  properties. 
Experimentally,  the  phenomena  of  flame  ignition,  propagation,  and  extinction  of  systematically 
chosen  fuel/oxidizer  mixtures  were  investigated.  The  studies  considered  both  gaseous  and  liquid 
hydrocarbons  in  view  of  their  importance  towards  the  accurate  description  of  the  combustion 
behavior  of  practical  fuels.  The  experimental  studies  expanded  notably  the  parameter  space  of 
existing  archival  fundamental  combustion  data.  The  detailed  modeling  of  the  experimental  data 
and  existing  ones  in  the  literature,  was  performed  by  using  well-established  numerical  codes  and 
reaction  and  transport  models  that  have  been  developed  based  on  a  novel  approach.  More 
specifically,  the  derivation  of  the  kinetics  models  was  based  on  the  following  four  specific 
objectives:  (a)  to  develop  and  test  a  suitable  array  of  mathematical  and  computational  tools  that 
can  be  used  to  quantify  the  Joint  rate-parameter  uncertainty  space  and  in  doing  so,  to  facilitate 
the  rational  design  of  reaction  models  suitable  for  any  given  target  fuel  or  mixture  of  these  fuels; 
(b)  to  examine  the  applicability  of  the  H2/CO/C,.4  reaction  model,  developed  under  the  current 
and  prior  AFOSR  support,  to  predict  the  phenomena  of  flame  ignition,  propagation,  and 
extinction  along  with  flame  structures,  (c)  to  examine  the  effect  of  inelastic  collision  on  binary 
diffusion  coefficients,  and  (d)  to  test  the  hypothesis  that  there  exists  a  critical  reaction  model  that 
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can  be  used  as  the  quantitative  kinetic  foundation  for  modeling  the  combustion  of  high  molecular 
weight  hydrocarbons.  The  approaches  to  meet  the  objectives  included  quantum  chemistry 
calculations,  reaction  rate  theory,  kinetic  modeling,  and  large-scale  optimization  using  the 
response  surface  method.  The  research  that  was  completed  not  only  produced  a  comprehensive 
reaction  model  suitable  for  air-breathing  propulsion  simulations,  but  it  also  generated  a  set  of 
new  approaches  and  computational  tools  for  rational  reaction  model  development  and 
optimization.  It  was  shown  systematically,  that  not  only  the  chemistiy  of  hydrocarbons  even  at 
the  Q  level  is  not  well  described  but  that  there  are  also  notable  uncertainties  associated  with  the 
Hj/Oj  chemistry  especially  when  attempts  are  made  to  predict  simultaneously  low -temperature 
and  high-temperature  flame  phenomena.  The  latter  is  a  rather  important  issue,  as  the  H^/O; 
kinetics  are  in  general  considered  to  be  rather  well  understood  and  are  a  very  important  subset  of 
the  kinetics  of  all  practical  hydrocarbon  fuels.  The  uncertainties  of  the  H2/O2  kinetics  can  have  a 
profound  effect  on  the  simultaneous  prediction  of  ignition  and  vigorous  burning  of  Hj  and 
hydrocarbons. 
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2.0  Experimental  Approach 

The  experiments  were  conducted  in  the  counterflow  configuration  as  shown  in  Rgs.  1  and  2, 
which  allows  for  the  establishment  of  steady,  laminar,  planar  premixed  or  non-premixed  flames 
and  the  measurement  of  laminar  flame  speeds,  as  well  as  ignition  and  extinction  limits.  From  a 
practical  point  of  view,  these  properties  are  measures  of  the  ability  of  flames  to  be  established 
and  sustained  under  the  extreme  conditions  encountered  in  engines.  From  a  fundamental  point 
of  view,  such  data,  being  sensitive  to  the  fuel’s  oxidation  chemistry,  can  be  used  to 
compile/validate/optimize  kinetic  models. 

An  array  of  counterflow  rigs  were  used  and  each  one  comprises  of; 

■  A  flow  system  that  is  appropriate  for  the  handling  of  gaseous  fuels  as  well  as  liquid 
fuels  in  their  vapor  phase 

•  Several  pairs  of  variable-diameter  heated  burners  opposing  each  other 

•  Either  a  digital  particle  Image  velocimetry  (DPIV)  or  a  laser  Doppler  velocimetry 
(LDV)  system  appropriate  for  velocity  measurements 

•  Several  thermocouple  system  appropriate  for  temperature  measurements. 

One  set  of  counterflow  burners  is  housed  into  a  stainless  steel  chamber,  shown  in  Fig.  3, 
which  can  operate  from  pressures  well  below  atmospheric  up  to  10  atm  and  is  equipped  with  the 
state-of-the-art,  remotely  controlled  functions. 

The  vaporization  system  consists  of  a  high-precision  syringe  pump  that  injects  the  fuels  into 
a  heated  chamber  surrounded  by  a  flow  of  heated  air  or  Nj.  The  flow  rates  of  all  gases  are 
measured  using  calibrated  sonic  orifices.  The  vaporization  system  has  undergone  several 
iterations  to  achieve  optimum  performance  and  reliability  for  heavy  liquid  fuels.  The 
vaporization  chamber  is  made  out  of  glass  chamber  to  prevent  fuel  cracking.  Additionally,  a 
quartz  nebulizer  with  a  flush  capillary  lapped  nozzle  was  integrated  into  the  system  in  order  to 
introduce  the  fuel  into  the  chamber  as  a  fine  aerosol;  the  details  are  shown  in  Fig.  4.  Thus,  rapid 
vaporization  is  achieved  at  relatively  low  chamber  temperatures,  allowing  steady-state 
vaporization  and  minimizing  the  chances  of  fuel  cracking.  To  prevent  fuel  condensation,  the  fuel 
delivery  line  is  wrapped  with  electrically  heated  tapes,  and  its  temperature  is  monitored  by 
properly  positioned  thermocouples.  The  burners  are  heated  with  ceramic  heating  jackets  to 
maintain  the  unburned  mixture  temperature,  Tu.as  high  as  500  K. 
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A  thermocouple  is  used  to  monitor  the  temperature  at  the  center  of  the  burner  exit.  Before 
each  experiment,  the  entire  system  is  allowed  to  equilibrate  and  reach  steady  state,  typically 
within  two  hours. 
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Figure  I  .  Flames  that  can  be  established  in  the  opposed-jet  configuration. 


Figure  2.  Schematic  of  a  typical  combustion  rig. 


6 


Figure  3.  High-pressure  chamber  equipped  with  counterflow  burners  and  a  state  of  the  art  DPIV  system. 


Figure  4,  Schematic  and  picture  of  the  updated  liquid  fuel  vaporization  system. 

Flame  ignition  is  achieved  by  heating  air  or  Nj  in  the  upper  burner  at  temperatures  well 
exceeding  1000  K.  To  achieve  this,  the  burner  is  replaced  by  a  quartz  tube  equipped  with 
independently  controlled  helicoidal  internal  (SiC)  and  external  (Kanthal  AI)  electric  heating 
elements  that  are  capable  of  raising  the  temperature  of  air  or  Nj  up  to  1600  K  (Fotache  et  al. 
1997).  A  temperature  control  ler/power  supply  is  used  to  operate  the  internal  heater.  Feedback  to 
the  temperature  controller  is  provided  by  a  B-type  thermocouple  (Pt-6%  Rh/Pt-30%  Rh) 
fX)sitioned  in  the  heated  air  or  Nj  stream  within  the  quartz  tube  close  to  the  surface  of  the  inner 
heater.  For  a  given  reactant  configuration  and  composition,  the  temperature  at  the  quartz  tube 
exit  is  increased  and  monitored  via  another  B-type  thermocouple,  until  flame  ignition  occurs,  at 
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which  point  the  tube  exit  temperature  is  defined  as  the  ignition  temperature,  T-,^.  Aspects  of  the 
ignition  system  are  shown  in  Fig.  5. 


Figure  5.  ignition  counterflow  burner  rig. 


The  measurements  included  axial  flow  velocities  and  temperature  measurements.  The  axial 
flow  velocities  were  measured  along  the  stagnation  streamline  using  either  LDV  or  DPIV.  The 
flow  was  seeded  with  0.3-p.m  diameter,  heated  silicon  oil  droplets  that  were  produced  by  a 
nebulizer.  A  typical  axial  velocity  profile  is  shown  in  Fig.  6.  The  local  strain  rate,  K,  is  defined 
as  the  maximum  absolute  value  of  the  velocity  gradient  Just  upstream  of  the  flame,  as  shown  in 


8 


Fig.  6  (Law  1988).  The  minimum  axial  flow  velocity  just  upstream  of  the  flame  is  defined  a 
reference  flame  speed,  as  shown  also  in  Fig.  6  (Law  1988). 


3.0  Numerical  Approach 

The  experiments  were  simulated  directly  by  including  detailed  descriptions  of  chemical 
kinetics  and  molecular  transport.  Laminar  flame  speeds,  5",  extinction  strain  rates,  K^^,  and 
were  computed  by  utilizing  CHEMKIN-based  codes  (Kee  at  al.  1985,  1986,  1988,  1989; 
Egolfopoulos  &  Campbell  1996).  Those  codes  have  been  upgraded  by  the  Pis  to: 

•  Account  for  any  type  of  combustion  mode  (Egolfopoulos  1994a;  Egolfopoulos  & 
Campbell  1996); 

•  Account  for  any  type  of  boundary  condition  (Egolfopoulos  1994a;  Egolfopoulos  & 
Campbell  1996,  Egolfopoulos  etal.  1997); 

•  Account  for  thermal  radiation  (Egolfopoulos  1994b;  Zhang  &  Egolfopoulos  2000); 

•  Solve  around  singular  turning-points  at  the  states  of  ignition  and/or  extinction 
(Egolfopoulos  &  Dimotakis  1998;  Dong  et  al.  2005;  Holley  et  al.  2006, 2009); 

•  Perform  mathematically  rigorous  sensitivity  analysis  of  all  dependent  parameters, 
including  S",  /sTo,,,  and  on  all  rate  constants  and  all  binary  diffusion  coefficients 
(Dong  et  al.  2005;  Holley  et  al.  2006, 2009). 

Additionally,  efficient  post-processing  codes  have  been  developed  to  perform  reaction  path 
analysis,  which,  when  combined  with  the  sensitivity  analyses,  provide  the  required  insight  into 
the  controlling  kinetic  pathways  for  each  phenomenon  that  is  investigated. 
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4.0  Results  and  Discussion 

4. 1  Development  of  a  new  non-linear  extrapolation  technique  for  the  determination  of  laminar 
flame  speeds 

To  determine  5°,  5u.ref  and  K  need  to  be  measured  first,  as  shown  in  Fig.  6.  Subsequently, 
plotting  Su,ref  against  K,  S°  could  be  determined,  in  principle,  by  extrapolating  linearly  the  Su.ref 
data  to  zero  stretch,  iC=  0  (e.g..  Law  1988).  However,  it  has  been  shown  (e.g.,  Tien  and  Matalon 
1991)  that  as  0  the  variation  of  5u.ref  may  not  be  linear  with  respect  to  K  because  of  the 
dilatation.  Therefore,  the  linearly  extrapolated  values  of  S°  have  been  shown  to  be  higher  than 
those  obtained  by  using  non-linear  extrapolation  (Davis  &  Law  1998)  or  by  other  more  direct 
methods  (Vagelopoulos  &  Egolfopoulos  1998). 

A  new  non-linear  extrapolation  technique  was  introduced.  Specifically,  the  variation  of  5u.ref 
with  K  is  computed  from  first  principles,  i.e.  using  detailed  description  of  chemical  kinetics  and 
molecular  transport,  over  a  range  of  K  for  which  there  is  no  heat  loss  to  the  burner  boundary.  In 
addition,  the  5“  value,  corresponding  to  K  =  0,  is  computed  in  order  to  anchor  the  5u,rervs.-/ir 
correlation.  The  computed  Su.ref  and  5°  values  are  then  fitted  as  a  polynomial  function  of  K.  It 
was  shown  that  small  but  finite  uncertainties  in  rate  constants  and  reactant  transport  coefficients 
have  minimum  effect  on  the  shape  of  the  5u,rcrvs.-A'  correlation.  Thus,  the  polynomial  function 
can  be  translated  vertically  to  fit  the  experimental  best  and  determine  5".  The  sensitivity  of 
the  shape  of  the  5'u,rervs.-/f  correlation  on  kinetics  was  tested  further  by  considering  various 
kinetic  models  that  predict  reasonably  well  S°,  and  it  was  determined  that  the  shape  of  the  non¬ 
linear  extrapolation  curves  were  affected  minimally  by  the  choice  of  the  model.  This  finding  is 
physically  reasonable.  More  specifically,  the  response  of  5ii.ref  to  fluid  dynamics  is  due  to  flame 
stretch  and  dilatation.  In  opposed-jet  flames  that  are  computed  using  similar  transport 
parameters  but  two  different  kinetic  models  that  predict  similar  S°  values,  the  effects  of  stretch 
and  dilatation  is  expected  to  be  similar.  Finally,  a  computed  5u.cBrvs.-^r  correlation  from  first- 
principle,  represents  the  experimental  data  better  compared  to  asymptotic  analysis  in  which  one- 
reactant,  one-step  kinetics  and  over-simplified  assumptions  for  the  transport  coefficients  are 
used. 

Similarly  to  previous  studies  invoking  non-linear  extrapolations  (e.g.,  Davis  &  Law  1998; 
Yang  et  al.  2009),  the  present  study  confirmed  that  linear  extrapolations  result  in  higher  5" 
values.  Figure  7  depicts  the  percentage  difference  between  linearly  and  non-linearly 
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extrapolated  laminar  flame  speeds  for  atmospheric  C5  -  C12  n-alkane  flames.  The  results 
illustrate  that  the  difference  is  about  5%  for  fuel-lean  mixtures  (equivalence  ratio  0<  I.O),  but  it 
rises  to  about  20%  for  fuel-rich  mixtures  {(p  >  1 .0).  Such  differences  are  rather  large  if  reaction 
kinetics  are  to  be  validated  against  S°  data.  The  increasingly  larger  discrepancies  towards  fuel- 
rich  mixtures  relates  to  their  sub-unity  Lewis  number  for  the  fuels  considered  herein,  given  the 
large  molecular  weight  discrepancy  between  the  fuel  and  oxygen  (e.g.,  Law  1988;  Tien  & 
Matalon  1991). 


Figure  7.  Average  percentage  difference  between  S° ’s  of  atmospheric  Cs-Ci^  n-alkane/air  mixtures  as  obtained 
by  linear  and  non-linear  extrapolations.  Non-linear  extrapolations  result  in  lower  S°  values. 


4.2  Molecular  transport  and  kinetics  effects  on  flame  propagation  and  extinction  in  liquid 
hydrocarbon  flames 

Lennard-Jones  12-6  potential  parameters  were  estimated  for  n-alkanes  and  7-aIkenes  with 
carbon  numbers  ranging  from  5  to  14.  These  parameters,  while  somewhat  uncertain,  form  the 
basis  for  a  comprehensive  analysis  of  the  responses  of  various  n-dodecane  flame  phenomena  to 
kinetic  and  transport  parameters.  The  logarithmic  sensitivity  coefficients  on  both  reaction  rates 
and  binary  diffusion  coefficients  were  calculated  for  laminar  flame  speeds  as  well  as  extinction 
strain  rates  for  premixed  and  nonpremixed  flames.  Results  revealed  that  die  response  of 
premixed  flames  is  generally  insensitive  to  the  fuel  diffusivity.  Coupled  with  the  lack  of 
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sensitivity  to  the  reaction  kinetics  specific  to  the  fuel  molecule,  the  computational  results  explain 
satisfactorily  the  low  sensitivity  of  the  phenomena  of  flame  propagation  and  extinction  to  the 
variation  of  the  size  of  n-alkane  fuels.  In  contrast,  the  extinction  strain  rates  of  nonpremixed 
flames  were  found  to  be  highly  sensitive  to  the  fuel  diflusivity,  in  that  larger  fuel  diffusivities 
render  the  flames  more  resistant  to  extinction.  Hence,  the  variation  of  of  nonpremixed 
flames  with  the  »-alkane  size,  as  observed  Holley  et  al.  (2007),  is  explained  satisfactorily  by  the 
variation  of  the  mass  diffusivity  of  the  fuel. 

Figure  8  illustrates  a  summary  of  the  logarithmic  sensitivity  coefficients  of  the  model 
responses  of  «-dodecane  flames  to  key  reaction  rate  parameters  and  binary  diffusion  coefficients. 
Laminar  flame  speeds  (5^-^)  and  extinction  strain  rates  for  premixed  flames  (iCexr^)  at  several 
representative  ^’s,  and  extinction  strain  rates  for  nonpremixed  flames  [iCexi-fuel+nitrogen  vs. 
oxidizer]  were  considered  for  conditions  that  are  specified  in  Table  1 . 

Several  observations  can  be  made  from  the  computational  results.  The  overall  influences  of 
reaction  rates  and  diffusion  fluxes  are  about  the  same  for  the  premixed  flame  responses.  The 
sensitivity  to  diffusion  originates  primarily  from  the  O2-N2  binary  diffusivity.  Kinetically,  the 
propagation  rates  of  the  premixed  flames  are  found  to  be  insensitive  to  the  chemistry  specific  to 
«-dodecane  as  shown  here,  since  its  decomposition  in  the  preheat  zone  of  the  flame  is  too  rapid 
to  be  rate-limiting  (You  et  al.  2009).  Rather,  the  sensitivity  spectrum  includes  only  the  rate 
coefficients  found  typically  in  H2/CO  and  small  hydrocarbon  combustion.  As  shown  in  Fig.  8, 
these  are  mainly  H  +  O2  O  +  OH,  CO  +  OH  ^  CO2  +  H,  and  H  +  O2  +M  — »  HO2  +M.  The 
fact  that  the  computed  responses  of  premixed  flames  are  insensitive  to  kinetic  or  diffusion 
parameters  specific  to  n-dodecane  suggests  that  the  measured  5"  and  ^ext  of  premixed  flames 
should  be  independent  of  the  nature  and  size  of  the  alkane  fuel  molecule.  Note  that  under  similar 
conditions,  the  decomposition  kinetics  of  higher  «-alkanes  is  similar  to  those  of  n-dodecane. 
Additionally,  for  smaller  n-alkanes  the  larger  diffusivities  would  cause  the  flame  propagation 
rate  to  depend  even  less  on  the  fuel  diffusion  rate. 
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Figure  8,  Logarithmic  sensitivity  coefficients  on  rate  constants  and  binary  diffysion  coefficients*  Light  bars: 
positive  values;  dark  bars:  negative  values. 


Flame 

Description 

Laminar  Flame  Saeed  (T..  =  403  K) 

s:  1.4 

^1.4 

5: 1.0 

0=1.0 

s:  0.7 

0=0.7 

Extinction  strain  rate  ofnremixed  flame 

Fuel/air  (^1.4)  vs*  N2 

Fuel/air  (^LO)  vs.  N2 

/r«.o.7 

Fuel/air  (^0.7)  vs.  N2 

Extinction  strain  rate  o  f  Jionpremixed  flame 

1-30  (air) 

Fuel/  N2=1/30  vs.  air 

1-15  (air) 

Fuel/  N2=1/15  vs.  air 

1-120(02) 

Fuel/  N2=1/120  vs.  O2 

K^,  1-60(0^ 

Fuel/  N2=1/60  vs.  O, 

Table  1,  Conditions  of  computed  results  shown  in  Figure  8.  For  all  caseSj  the  fuel  is  w-dodecane,  the  pressure  is 
atmospheric,  and  the  temperature  of  the  fuel-containing  stream  is  403  K*  The  oxidizer  is  air  for  premixed  flames, 
and  air  or  O2  for  nonpremixed  flames.  For  nonpremixed  flames  the  fuel  stream  consists  of  rt-dodecane/nitrogcn 
mixtures  with  varying  molar  ratios* 
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4.3  Quantification  of  the  effect  of fuel  blending  on  the  flame  response 

iCext’s,  and  fjgn’s  were  computed  Laminar  flame  speeds  as  well  as  ignition  and 
extinction  limits  of  mixtures  of  air  with  hydrogen,  carbon  monoxide,  and  C1-C4  saturated 
hydrocarbons  were  studied  both  experimentally  and  numerically.  The  fuel  mixtures  were  chosen 
in  order  to  gain  insight  into  potential  kinetic  couplings  during  the  oxidation  of  fuel  mixtures. 
The  experiments  were  modeled  using  the  USC  Mech  II  kinetic  model  (Wang  et  al.  2007).  It  was 
determined  that  when  hydrocarbons  are  added  to  hydrogen  flames  as  additives,  flame  ignition, 
propagation,  and  extinction  are  affected  in  a  counterintuitive  manner.  More  specifically,  it  was 
found  that  by  substituting  methane  by  propane  or  n-butane  in  hydrogen  flames,  the  reactivity  of 
the  mixture  is  reduced  both  under  pre-ignition  and  vigorous  burning  conditions.  This  behavior 
stems  from  the  fact  that  propane  and  n-butane  produce  higher  amounts  of  methyl  radicals  that 
can  readily  recombine  with  atomic  hydrogen  and  reduce,  thus,  the  rate  of  the  H  +  O2  ^  O  +  OH 
branching  reaction.  The  kinetic  model  predicts  closely  the  experimental  data  for  flame 
propagation  and  extinction  for  various  fuel  mixtures  and  pressures.  On  the  other  hand,  it 
underpredicts,  in  general,  Tign’s. 

Additionally,  5"’s  of  mixtures  of  air  with  80%  n-dodecane  +  20%  methylcyclohexane  and 
80%  n-dodecane  +  20%  toluene,  on  a  per  volume  basis,  were  determined  over  a  wide  range  of 
0’s,  at  atmospheric  pressure,  and  ~  403  K.  The  choice  of  the  fuel  blends  was  dictated  by  their 
anticipated  compositions  in  surrogates  of  jet  fuels.  Phenomenological  analysis  revealed  that 
5“’s  of  binary  fuels  mixtures  can  be  estimated  using  5"’s  and  adiabatic  flame  temperatures  of 
flames  of  the  neat  components.  The  propagation  rates  of  flames  of  the  various  binary  fuels 
blends  were  computed  using  detailed  descriptions  of  chemical  kinetics  and  molecular  transport 
and  were  found  to  be  in  good  agreement  with  the  estimations.  Although  the  initial  fuel  initial 
consumption  pathways  and  the  resulting  intermediates  and  radicals  may  be  different  for  each 
neat  component,  the  propagation  of  flames  of  binary  fuels  was  found  to  be  mostly  sensitive  to 
the  flame  temperature  through  its  influence  on  the  main  branching  reaction  H  +  O2  ^  OH  +  O. 
Thus,  kinetic  couplings  resulting  from  the  presence  of  two  different  fuels  appear  to  have  minor 
effect  on  flame  propagation. 
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4.4  Derivation  of  an  extensive  experimental  flame  database  for  C5-C12  n-alkane,  cyclohexane, 
and  mono-alkylated  cyclohexane  flames 

5°’s  and  of  premixed  C5-C12  «-alkane  flames  were  determined  over  a  wide  range  of  0’s, 
at  atmospheric  pressure,  and  elevated  Tu’s.  5"’s  were  obtained  using  the  aforementioned  in 
Section  4.1  non-linear  extrapolation  technique.  Compared  to  linearly  extrapolated  values,  the 
laminar  flame  speeds  obtained  using  non-linear  extrapolations  were  found  to  be  1  to  4  cm/s 
lower  depending  on  <p.  It  was  determined  that  S"’s  of  all  /i-alkane/air  mixtures  considered  in  this 
investigation  are  similar  to  each  other  and  sensitive  largely  to  the  H2/CO  and  C1-C4  hydrocarbon 
kinetics.  Additionally,  the  resistance  to  extinction  decreases  as  the  fuel  molecular  weight 
increases.  Simulations  of  the  experiments  were  performed  using  the  recently  developed  JetSurF 
0.2  reaction  model  (Sirjean  et  al.  2008)  consisting  of  194  species  and  1459  reactions.  The 
experimental  5"’s  were  predicted  with  good  accuracy  for  all  the  «-alkanes.  The  experimental 
text’s  were  predicted  closely  by  the  model  for  fuel-lean  mixtures.  For  stoichiometric  and  fuel- 
rich  mixtures,  the  computed  Ke^s  were  approximately  10%  lower  than  the  experimental  values. 
Insights  into  the  physical  and  chemical  processes  that  control  the  response  of  n-alkane  flames 
were  provided  through  detailed  sensitivity  analyses  on  both  reaction  rates  and  binary  diffusion 
coefficients.  Figures  9  and  10  depict  the  experimental  and  computed  5°’s. 

5"’s  were  measured  in  the  counterflow  configuration  for  cyclohexane/air, 
methylcyclohexane/air,  ethylcyclohexane/air,  n-propylcyclohexane/air,  and  n- 
butylcyclohexane/air  mixtures,  at  atmospheric  pressure,  Tu  =  353  K,  and  for  a  wide  range  of  0’s. 
It  was  determined  that  cyclohexane/air  flames  propagate  somewhat  faster  than  mono-alkylated 
cyclohexane/air  flames.  Flames  of  mono-alkylated  cyclohexane  compounds,  i.e.  methyl-,  ethyl-, 
n-propyl-,  and  «-butyl-cyclohexane,  were  found  to  have  similar  5"’s,  suggesting  that  the 
different  alkyl  groups  have  a  secondary  effect  on  flame  propagation.  The  experimental  data 
were  modeled  using  JetSurF  1.1  (Sirjean  et  al.  2009),  a  detailed  reaction  model  for  the 
combustion  of  cyclohexane  and  its  derivatives,  and  satisfactory  agreements  were  found.  Based 
on  the  analysis  of  the  modeling  results,  the  somewhat  lower  rates  of  mono-alkylated  cyclohexane 
flame  propagation  were  attributed  to  the  greater  production  of  propene  and  allyl  and  the 
increased  H-atom  scavenging  by  these  Cj  intermediates.  Though  these  fuel-specific  reaction 
kinetic  features  do  not  limit  the  overall  oxidation  rates,  the  distribution  of  the  cracked  products 


Laminar  Flame  Speed,  cm/s  Laminar  Flame  Speed,  cm/s 
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exert  influences  on  flame  propagationj  leading  to  the  subtle  differences  in  5^’s  that  were 
observed  experimentally.  Figure  1 1  depicts  the  experimental  and  computed 


Figure  9,  Experimentally  and  computed  at  =  353  K,  (a)  «-C5Hi2,  (b)  (c)  n-C7Hi6,  and  (d) 

w-CbH I g/air  mixtures.  Symbols:  present  experimental  data  obtained  with  non-linear  extrapolation;  Lines: 
simulation  results  using  the  JetSurF  0.2  kinetic  model. 
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Figure  10*  Experimental  !y  and  computed  at  =  403  K.  (a)  w-C9H2f)/air  mixture,  (•)  Present 
experimental  data  obtained  with  non-linear  extrapolation;  (b)  n-CioH22/air  mixture,  (•)  Present  experimental 
data  obtained  with  non-linear  extrapolation;  (o)  Experimental  data  of  Kumar  &  Sung  (2007);  Lines: 
simulation  results  using  the  JetSurF  0,2  kinetic  modeL  The  error  bars  indicate  2-fj standard  deviations. 
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Figure  IL  Experimental  and  computed  of  methyl-,  ethyl,  /t-propyU  and  rt-butyl-cyclohexane/air  flames  at 
p  =  I  atm.  Symbols:  (H)  present  experimental  data  at  7'^  =  353  K,  (A)  Dubois  et  aL  (2009),  Lines:  ( — )  simulation 
using  Model  1  (JetSurF  1 .1)  at  7^  =  353  K,  ( — )  simulation  using  Model  I  (JetSurF  1 . 1)  at  7u  =  403  K.  The  error  bars 
indicate  2- a  standard  deviations. 


17 


4.5  Fundamental  kinetic  studies  of  elementary  reactions  critical  to  a  predictive  combustion 
reaction  model 

The  kinetics  of  the  reaction 

CO  +  HO2*  ->  CO2  +  'OH  (R1 ) 

was  studied  using  a  combination  of  ab  initio  electronic  structure  theory,  transition  state  theory, 
and  master  equation  modelling.  The  potential  energy  surface  was  examined  with  the  coupled 
cluster  CCSD(T)  method  and  complete  active  space  CASPT2  methods.  The  classical  energy 
barriers  were  computed,  as  shown  in  Fig.  12.  The  energy  barriers  were  found  to  be  18  and 
19  kcal/mol  for  CO  +  HOj*  addition  following  the  trans  and  cis  paths,  respectively. 


rs*35.1  rS7  35.7 


Figure  12.  Potential  energy  diagram  for  CO  +  H02*  — *■  products.  For  CO  +  H02*  “*-002  +  'OH,  the  energy 
values  were  determined  using  the  CCSD(T)/CBS  method,  and  included  zero-poini  energy  corrections.  For 
CO  +  H02*  ^  HC*0  +  O27  the  energy  values  were  taken  from  Martinez-Avila  et  al  (2003)* 

For  the  cis  path,  rate  constant  calculations  were  carried  out  with  canonical  transition  state  theory. 
For  the  trans  path,  master  equation  modelling  was  also  employed  to  examine  the  pressure 
dependence.  Special  attention  was  paid  to  the  hindered  internal  rotations  of  the  H00C*0  adduct 
and  transition  states.  The  theoretical  analysis  showed  that  the  overall  rate  coefficient  is 
independent  of  pressure  up  to  500  atm.  The  computed  rate  coefficient  was  found  to  be 
substantially  smaller  than  some  of  the  rate  recommendations  in  the  literature,  as  shown  In  Fig. 
13.  Based  on  the  current  theoretical  calculations,  the  following  rate  expression  was 
recommended  for  reaction  R1 

k^iem^  I  mol  -s)  =  \  .57jcl0^r^’‘ 
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for  300  <  rs  2500  K  with  the  uncertainty  factor  equal  to  8,  2  and  1.7  at  temperatures  of  300, 
1000,  and  2000  K.,  respectively. 
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Figure  1 3.  Rate  coefficient  of  reaction  RK 

The  gas-phase  reaction  of  benzene  with  0(^P)  is  of  considerable  interest  for  modeling  of 
aromatic  oxidation,  and  also  because  there  exist  fundamental  questions  eonceming  the 
prominence  of  intersystem  crossing  in  the  reaction.  While  its  overall  rate  constant  has  been 
studied  extensively,  still  there  are  significant  uncertainties  in  the  product  distribution.  The 
reaction  proceeds  mainly  through  the  addition  of  the  O  atom  to  benzene,  forming  an  initial  triplet 
diradical  adduct,  which  can  either  dissociate  to  form  the  phenoxy  radical  and  H  atom,  or  undergo 
intersystem  crossing  onto  a  singlet  surface,  followed  by  a  multiplicity  of  internal  isomerizations, 
leading  to  several  possible  reaction  products.  In  this  work,  the  product  branching  ratios  of  the 
reaction  was  studied  over  the  temperature  range  of  300  to  1000  K  and  pressure  range  of  1  to 
lOTorr.  The  reactions  were  initiated  by  pulsed-laser  photolysis  of  NO2  in  the  presence  of 
benzene  and  helium  buffer  in  a  slow-flow  reactor,  and  reaction  products  were  identified  by  using 
the  multiplexed  chemical  kinetics  photoionization  mass  spectrometer  operating  at  the  Advanced 
Light  Source  (ALS)  of  Lawrence  Berkeley  National  Laboratory.  Phenol  and  phenoxy  radical 
were  detected  and  quantified.  Cyclopentadiene  and  cyclopentadienyl  radical  were  identified 
directly  for  the  first  time.  Finally,  ab  initio  calculations  and  master  equation/RRKM  modeling 
were  used  to  reproduce  the  experimental  branching  ratios,  yielding  pressure-dependent  rate 
expressions  for  the  reaction  channels,  including  phenoxy  +  H,  phenol,  and 
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cyclopentadiene  +  CO,  as  shown  in  Fig.  14.  These  computed  rate  coefficients  were  used  to 
update  the  USC_Mech  II  reaction  model  (Wang  et  al.,  2007)  so  that  predictions  of  aromatics 
combustion  are  improved. 


laooK/r  1000K/T  ioook/t 

Figure  14.  Arrhenius  plots  for  the  reaction  channels  leading  to  (a)  phenol,  (b)  phenoxy  +  H,  (c)  cyclopentadiene  + 
CO,  (d)  benzene  oxide,  (e)  2,4-cyclohexadienone,  and  (f)  2,5-cyclohexadienone.  Symbols  are  rates  computed  with 
Monte  Carlo  RRKM/master  equation  modeling  and  lines  represent  fitted  Arrhenius  expressions. 
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4.6  Refinement  ofUSC  Mech  II  to  predict  H^CO  combustion  at  pressures  up  to  500  atm 

The  high-pressure  oxidation  of  inert-diluted  CO/O2  mixtures  doped  with  150-200  ppm  of  Ha 
was  studied  behind  reflected  shock  waves  in  the  University  of  Illinois  at  Chicago  high-pressure 
single  pulse  shock  tube.  The  experiments  were  performed  for  stoichiometric  (0=1)  and  fuel 
lean  (0=0.5),  over  the  temperature  and  pressure  ranges  of  1000-1500  K  and  21-500  bar 
respectively.  Stable  species  sampled  from  the  shock  tube  were  analyzed  by  gas  chromatographic 
(GC)  and  GC/mass  spectrometric  techniques.  The  experimental  data  were  simulated  using  the 
use  Mech  II  kinetic  model  (Wang  et  al.  2007).  These  simulations  showed  that  within  the 
experimental  uncertainty,  USC  Mech  11  was  able  to  capture  the  experimental  trends  for  the  lower 
pressure  data  sets  (average  nominal  pressures  of  24  and  43  bars).  However  the  model 
underpredicted  the  CO  and  O2  decay  and  subsequent  CO2  formation  for  the  higher-pressure  data 
sets  (average  nominal  pressures  of  256  and  450  bars).  The  elevated  pressure  data  sets  span  a 
previously  unmapped  regime  and  have  served  to  probe  HO2  radical  reactions,  which  appear  to  be 
among  the  most  critical  reactions  under  the  experimental  conditions  that  were  considered.  With 
updated  rate  parameters  for  a  key  HO2  radical  reaction,  USC  Mech  II  was  shown  to  reconcile  the 
elevated  pressure  data  sets  thereby  extending  the  model  predictability  to  an  extreme  range  of 
conditions. 
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Figure  15,  Experimental  data  for  CO  {•)  and  CO^  («)  measured  by  shock  heating  of  a  mixture  containing  500  ppm 
CO,  200  ppm  H2,  O2  and  argon  with  the  equivalence  ratio  equal  to  t,0.  The  computational  results  are  □:  base  USC 
Mech  II,  x:  updated  USC  Mech  IL 
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4.7  Development  of  detailed  and  simplified  models  for  the  oxidation  of  CyCn  n-alkanes 

A  detailed  kinetic  model  was  proposed  for  the  combustion  of  n-alkanes  up  to  n-dodecane  and 
for  temperatures  above  850  K.  The  model  was  validated  against  experimental  data,  including 
fuel  pyrolysis  in  plug  flow  and  jet-stirred  reactors,  laminar  flame  speeds,  and  ignition  delay  times 
behind  reflected  shock  waves,  with  n-dodecane  being  the  emphasis.  Selected  model  validation 
results  are  shown  in  Figs.  16  and  17.  The  model  was  used  as  the  basis  for  JetSurF  0.1  and  1.0 
with  a  wider  range  of  validation  tests  reported  on  the  JetSurF  web  site  (Sirjean  et  al.  2008;  2009) 
and  in  Figs.  9  and  10. 
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Figure  16.  Experimental  (symbols:  Herbmet  et  al.  2007) 

and  computed  (solid  line:  detailed  model;  dashed  line: 
simplified  model)  conversions  of  /i-Ci2H26  as  a  function 
of  the  residence  time  for  the  pyrolysis  of  2%  71-C12H26  in 
He  in  ajei-stirred  reactor  (p  -  1  atm). 


Figure  17.  Experimental  (symbols)  and  computed  (solid 
line:  detailed  model;  dashed  line:  simplified  model  not 
shown  here)  ignition  delay  times  for  «-CioH22  oxidation 
behind  reflected  shock  waves. 


Analysis  of  the  computational  results  revealed  that  for  a  wide  range  of  combustion 
conditions,  the  kinetics  of  fuel  cracking  to  form  smaller  molecular  fragments  is  fast  and  may  be 
decoupled  from  the  oxidation  kinetics  of  the  fragments.  Subsequently,  a  simplified  model 
containing  a  minimal  set  of  4  species  and  20  reaction  steps  was  developed  to  predict  the  fuel 
pyrolysis  rate  and  product  distribution.  Combined  with  the  base  C1-C4  model,  the  simplified 
model  predicts  fuel  pyrolysis  rate  and  product  distribution,  laminar  flame  speeds,  and  ignition 
delays  as  close  as  the  detailed  reaction  model. 
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4.8  Development  of  the  method  of  uncertainty  minimization  by  polynomial  chaos  expansion 
The  ability  of  a  reaction  model  to  predict  the  combustion  behavior  of  a  fuel  relies  on  the 
rigorous  quantification  of  the  kinetic  rate  parameter  uncertainty.  Although  the  accuracy  of  a 
detailed  kinetic  model  may  be  ensured,  in  principle,  by  a  multi-parameter  optimization,  the 
inherent  uncertainties  in  the  fundamental  combustion  targets  used  for  optimization  cause  the 
derived  optimized  model  to  be  characterized  by  a  finite  kinetic  parameter  space.  In  this  work, 
spectral  expansion  techniques  were  developed  and  employed  to  quantify  these  uncertainties.  The 
resulting  method,  called  the  Method  of  Uncertainty  Minimization  by  Polynomial  Chaos 
Expansion  (MUM-PCE)  was  used  to  analyze  a  detailed,  H2/CO/C1-C4  kinetic  model  for  ethylene 
combustion.  Uncertainty  was  quantified  for  both  the  as-compiled  model  and  the  optimized 
model,  and  was  propagated  subsequently  into  a  wide  variety  of  combustion  experimental 
conditions.  Examples  for  laminar  flame  speeds  and  shock  tube  ignition  delays  are  shown  in 
Figs.  18  and  19. 
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Figure  18.  Experimental  data  (o:  Egolfopoulos  et  al. 
(1990),  Jomaas  ei  al,  (2005),  0:  Hassan  et  al.  {1998)) 
and  2cr  uncertainty  bands  computed  for  the  laminar 
flame  speed  of  ethylene-air  mixtures  aXp  =  \  atm. 
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Figure  19-  Experimental  data  (symbols,  Brown  and 
Thomas,  1999)  and  computed  (dashed  line:  as- compiled, 
unoptimized  model;  solid  line:  optimized  model,  bands:  2o 
model  uncertainties)  ignition  delay  lime  of  1  %C2H4-3%02- 
Ar  mixtures  behind  reflected  shock  waves. 
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4.9  Gas-kinetic  theory  analysis  of  scattering  between  gas  and  small  particles  with  application  in 
the  transport  properties  of  high  molecular  weight  hydrocarbons  in  dilute  gases 

A  theoretical  study  was  conducted  to  examine  the  uncertainties  in  the  transport  theory  of 
molecular  clusters  and  small  particles  in  laminar  flow  regime  (Li  and  Wang  2003a,  2003b, 
2005).  The  theory  features  a  rigorous  gas-kinetic  theory  analysis,  considers  the  effect  of  non- 
rigid  body  collision,  and  has  been  shown  to  reproduce  the  Chapman-Enskog  theory  of  molecular 
transport  in  the  small  particle  size  limit,  Epstein's  model  of  particle  drag  in  the  rigid-body  limit, 
and  the  Stokes-Cunningham  equation  for  the  drag  of  micrometer  size  particles,  A  key 
uncertainty  in  the  application  of  the  generalized  theory  lies  in  the  uncertainty  of  momentum 
accommodation  of  molecules  or  particles  in  the  size  range  of  a  few  nanometers.  The  cause  for 
this  uncertainty  was  attributed  to  a  transition  in  the  dominant  mode  of  gas-particle  interactions 
from  specular-like  scattering  to  diffuse  scattering.  As  an  example.  Fig.  20  depicts  this  transition 
by  comparing  the  theoretical  reduced  collision  integrals  for  polyethylene  glycol  (PEG)-air 
molecule  interactions  and  the  experimental  ones  derived  by  Fernandez  de  la  Mora  and  coworkers 
(Nasibulin  et  al.,  2002;  Ude  and  Fernandez  de  la  Mora,  2003).  It  was  concluded  that  in  order  to 
apply  accurately  a  generalized  transport  theory  for  large  molecules,  a  substantially  larger  set  of 
experimental  mobility  and  difflisivity  has  to  be  available. 


R{nm) 

Figure  20.  Variations  of  the  reduced  collision  integrals  as  a  function  of  the  radius  of  PEG  nanoparticles.  Symbols 
are  experimental  data  observed  in  air  at  room  temperature  (diamonds:  Nasibulin  et  al.,  2002;  circles:  Ude  and 
Femdndez  de  la  Mora,  2003). 
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Technical  Interactions  with  AFRL  Researchers 

Several  meetings  with  Dr.  Tim  Edwards  of  AFRL  took  place  during  the  reporting  period. 
The  purpose  of  these  meetings  was  to  discuss  the  measurements  on  the  propagation,  ignition,  and 
extinction  limits  of  flames  of  neat  hydrocarbon  fuels  as  well  as  a  variety  of  petroleum -derived 
and  synthetic  jet  fuels  that  were  supplied  to  the  USC  Combustion  and  Fuels  Laboratory  by  Dr. 
Edwards.  The  discussions  has  focused  on  both  the  experimental  and  modeling  work  pertaining 
to  jet  fuels  and  their  surrogates,  as  well  as  ongoing  work  on  smaller  hydrocarbons  in  the  C1-C4 
range. 


Other  Technical  Interactions 

A  close  interaction  was  established  with  Professors  Ron  Hanson,  Tom  Bowman  and  Heinz 
Pitsch  of  Stanford  University,  Professor  C.K.  Law  of  Princeton  University,  Professors  Nick 
Cemansky  and  David  Miller  of  Drexel  University,  Professor  Angela  Violi  of  University  of 
Michigein,  Dr.  Wing  Tsang  of  National  Institute  of  Standards  and  Technology,  and  Dr.  Meredith 
Colket  of  United  Technologies  Research  Center  on  the  experimental  and  modeling  aspects  of 
this  research.  Additionally,  collaborative  research  efforts  were  initiated  with  Dr.  Stephen 
Klippenstein  of  the  Argonne  National  Laboratory  and  Dr.  Craig  Taatjes  of  the  Sandia  National 
Laboratories 


30 


Inveniions,  Technology  Transitions  and  Transfers 

Part  of  the  work  described  in  this  report  resulted  in  the  following  technology  transition  and 
transfer: 


Performer:  Dr.  Hai  Wang,  University  of  Southern  California,  (213)  740-0499 

Customer:  Dr.  Meredith  Colket,  United  Technology  Research  Center,  41 1  Silver  Lane 

East  Hartford,  CT  061  OS 
Phone:  (860)  610-7000 

Result:  Reaction  model  for  predicting  aromatics  and  soot  formation  in  fuel  rich  flames 

Application:  Low-emission  gas  turbine  engines 


